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Abstract 

The leading nonperturbative QCD corrections to the one gluon exchange quark-quark, 
quark-antiquark and qq pair-excitation potentials are derived by using a covariant form of 
nonlocal two-quark and two-gluon vacuum expectation values. Our numerical calculation 
indicates that the correction of quark and gluon condensates to the quark-antiquark potential 
improves the heavy quarkonium spectra to some degree. 
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1. Introduction 

A remarkable success of the potential model in explaining hadronic spectra and hadronic 
properties has been made ]]], 0, Inspired by QCD, various potential forms 0], in which 
both the Coulomb potential arising from the one gluon exchange (OGE) mechanism and 
confinement potential are employed, can give results which are reasonably consistent with 
the data at the charm and bottom energy scales within a certain error. 

Along a parallel direction to the QCD sum rule approach ||, Shen et al. || derived the 
leading nonperturbative QCD correction to the OGE quark-quark potential and investigated 
possible effects of the nonperturbative contribution in heavy quarkonium spectra. In this 
paper, we will extend the nonperturbative calculation of Ref. || to the nonperturbative 
correction to the quark-antiquark and qq pair-excitation potentials in a covariant form of the 
nonlocal two-gluon vacuum expectation value (VEV). 

According to some new features which appear in the newly obtained quark-quark potential 
including the nonperturbative correction, there are some intended applications of the mod- 
ified OGE potentials. Firstly, the modified quark-quark potential can be used to study the 
nonperturbative effect in the spectra of J/\I/ and T, especially to improve the spin splitting 
for these system. Because the (qq) excitations give significant contributions to the nucleon 
and delta masses, we also hope that modified qq excitation potential can be employed to 
investigate possible nonperturbative effects in these hadronic properties. In addition, the 
nonperturbative effect in decay constants of the pseudoscalar mesons can also be consid- 
ered by involving the vacuum condensate effect. Even some long-standing problems in light 
baryon spectroscopy such as the energy level order between the positive- and negative- parity 
partner states, in particular, Roper resonance puzzle, and the baryon spin-orbit structure 
puzzle-the absence of spin-orbit effects in baryon spectrum can partly be understood along 
this direction. 

In the nonrelativistic reduction of the perturbative OGE potential, if the CVC theorem 
is respected, i.e. q tJ "u q (pi)^f ll u q (p' 1 ) = with q = pi — p[, the term q^/q 2 in the S-matrix 
for a quark-quark scattering vanishes, and the quark-quark interaction is not affected by the 
gauge dependence of the gluon propagator. 

However, the gauge issues exist when both the gluon propagators and the external QCD 
vacuum gluon fields appear in the S-matrix for the nonperturbative quark-quark scattering. 
For the gluon propagators, we will adopt Landau gauge as done in Refs. Q and since it 
has been pointed out in Ref. [|| that a Landau gauge calculation for an on-shell amplitude 
can yield a correct gauge-invariant result; For the external QCD vacuum gluon fields, the 
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fixed-point gauge |J is generally adopted, i.e., 



x^(x) = 0. (1) 
Hence the nonlocal two-gluon VEV is || |l^] 



(0\B a Jx)B b v (y)\0) = ^V(0|G" G^|0) 



48(N2-l) xPya{gpa9flu ~ 9p^)(0\G 2 \0) + • • • , (2) 



where 



<0|C|0> = <0|G^Gf|0). (3) 

The expansion of (§) obviously violates the translational invariance since the right hand 
side(RHS) of (||) is a function of xy instead of (x — y). 



In Ref. [11], we pointed out that there is an uncertainty in the evaluation of the gluon 
condensate contribution to a three-point nonperturbative vertex provided that the fixed-point 
gauge of the QCD vacuum gluon fields is adopted. We found in Ref. |l2| that, by using a 
covariant form of the nonlocal two-gluon VEV, the gluon condensate correction to a fermionic 
three-point vertex can be unambiguously deduced. 

In order to avoid the ambiguity in the fixed-point gauge, it is better to derive the conden- 
sate correction to the OGE potentials by means of the covariant form of (0\B° l (x)B b (y)\0) 



given by Bagan et al. p3| , 

A_, r „ 9 1 

+ , (4) 



(0\B^x)Bl(y)\0) = -J^rrf 



2 2 

{x - y) - -{x - - y) L 



where 

a. 5(0|G 2 |0) 

b ~ 288 ' (5) 
Deriving the condensate corrections to the perturbative OGE potentials by means of the 
covariant form of nonlocal two-quark and two-gluon VEVs is the main purpose of this paper. 

In the following three sections, using the covariant form of nonlocal two-gluon VEV [O ], 
we derive the leading nonperturbative QCD corrections to the perturbative OGE quark- 
quark, quark-antiquark and qq pair-excitation potentials, respectively. Finally, we present a 
brief discussion and summary with a preliminary numerical result. 

2. Quark-Quark Interaction with Correction from Nonperturbative QCD 
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The nonrelativistic reduction method is used to extract the nonperturbative correction to 
the OGE potentials. The Feynman diagrams that will be used to evaluate the nonperturbative 
correction to quark-quark interaction are shown in Fig. 1. For the scattering of the two 
quarks of different flavors, by using the covariant form of two-gluon VEV of Ref. ||13[| , the 
contribution of the Feynman diagram 1(a) to the S-matrix can be written as 



x / ctk 



TV? 



i \ 9pa9lm 5 9 P l9 ° m j dhdk m 



o 2 



5\k) 



x<W/«'cd (2q + k)Pg*' x + (-q - 2kf g px + (k 



)V M ' 



x/, 



c'b'd' 



(2q + kYg x ' u ' + (-q + kYg v '° + {-2k 



A' v'a 



9*° D b X(q)- 



(6) 



Eq. (pt) can be rewritten as a relation between the corresponding quark potential: 



768ira s N c S .qgep 



(?)> 



(7) 



l)g 4 qq 

where V® GEP (q) is the usual perturbative quark-quark potential arising from the OGE mech- 
anism. Performing Fourier transformation to Vi^(q) of (^), we obtain the contribution of 
the Feynman diagram 1(a) to the quark-quark potential in the coordinate representation, 



A-^xf + A-l\x\ 



(8) 



where 



A, 



32N a S 



1) y raim 2 J 



(9) 



A x 



48N C S \ {m l + m 2 ) 2 1 
o^iV^ — \)m\mi \ m\m2 2 



1 

+ 2 



(2 + ^ + (2 + ^ 
mi m2 



(x x p) 



(10) 



with n = x/\x\. Fig. 1(b) makes no contribution, i.e., Uu^^x) = 0, because the vacuum 
cannot possess any quantum number such as color. 
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By means of the two-quark VEV |]1 
the S-matrix can be expressed as 



the contribution of the Feynman diagram 1(c) to 



X J d A k5 4 (k)(0\q f q f \0) 



-ig^-(p'2)i v ^ + {P2)W{A)i^ + (Pi)\ 



1 



f 



d 



r 



4N C 16N C 1 dk 



Y-Ls(q + k) 
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2 w 

where the next-to-leading-order term in the full coefficient of the (qq) component of the non- 
perturbative two-quark VEV [|l(]] is retained. There exists a similar expression for Fig. 1(d). 
It is easy to find that the contributions of the Feynman diagrams 1(c) and 1(d) to the effective 
potential are the same. Therefore, we can extract the effective potential for the scattering of 
two quarks of different flavors as shown in Fig. 1(c) and Fig. 1(d) in the coordinate represen- 
tation 



where 



Ui( c )(x) = ^i(d)(^) 



8(t) 



A 1 A 2 2 

— ; — Tra: 



Ci\x\ + C_i|f| _i + C- 3 \x 



(/)|^f|-l N \ „-m f \x\ 



c)!>\x\ 



(12) 



\ m\mi J JS c m,f 



(13) 



C-i 



2 y <0|g/g/|0> 

y m f 
i 



N c 



— — (<7i • n)(<? 2 -n) + - (— + — 
lm\mi 8 \mi m,2 



+ 



2mirri2 



(2 + ^-)a 1+ (2 + ^-)a 2 ]-(xxp)\, (14) 



3 ^ ((%/<j/|0) 
mim2-/V c rrij 



(15) 



???. 



/ 



m\m2 



1/1 1 

8 \77tl 1712 



+ 



1 



4mim2 



[(<?i • <x 2 ) - (<?i • n)(a 2 -n)}}, 



(16) 
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and 



1 <0|g /g/ |0) 



f 



-r— [(^1 • ^2) + (<?i • n){a 2 ■ ft)] + j f — + — ) 

4mim2 4 \mi ^2/ 



1 



2 + — )<?!+ (2 + — )«T 2 



(fxjj) 



(17) 



2mi77l2 

The total quark-quark effective potential is finally obtained by summing up the contributions 
of all the corresponding diagrams including perturbative and nonperturbative ones: 



U q Jx) = U° G ™(x) + U™(x) 



(18) 



where C/^ p (x), the nonperturbative correction to the perturbative quark-quark potential due 
to the quark and gluon condensates, can be expressed as 



^1(8)0) + U 1{h) {x) + U 1{c) (x) + U 1(d) (x 
S(t) 



A 1 A 2 2 

— ; 7TOU 



+2C_ 3 |fr 3 + 2 



[A 3 \x\ 3 + (A x + 2Ci) |x| + 2C_i|fr 1 

e -m f \x\ 



(19) 



Formally, Eq. fllq ) holds not only for the </(/-, but also for the qq- and qq- interactions. Note, 
however, that the color generators for an antiquark are given by — X T , i.e., 



U^ ect (x) = M*)Ia ? ash-a ?( a«F ( 2 °) 

and 

Uqq(x) = C/'gg(x)| A a A a^ {A a ) T {A a ) T. (21) 

For an interaction between a quark and an antiquark, if the quark and antiquark are of 
the same kind of quark fields, then not only the direct scattering but also the annihilation 
mechanism should be taken into account. The detail discussion is given in the following 
section. 

3. Quark- Antiquark Annihilation with Correction from Nonperturbative QCD 

In the above section, the direct interaction potential between the quark and antiquark was 
given. When the quark and antiquark have the same kind of flavor, then the annihilation of a 
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quark and an antiquark is possible. The annihilation diagrams including the nonperturbative 
correction shown in Fig. 2 should be further considered. We can obtain the contributions of 
these Feynman diagrams to the quark-antiquark annihilation potential by means of the same 
procedure as above. Taking Fig. 2(a) as an example, the S-matrix for this diagram can be 
obtained by making the substitutions 



(22) 



in Eq. (^). In the calculation of color part and spin part in the S-matrix for diagrams in 
Fig. 2, it is valuable to note that 



N*-l 



Tf ( A )«'«l A W ^ ^ d a'a0(3f3> 



( A<1 )/3a(A a )a' 



iV 



AT2 



(23) 



a=l a=l 

for SU(iV) generators X a (a = 1, 2, • • • , iV 2 — 1), with 7V=3. 2 for the color and spin generators, 
respectively. Thus, the effective potential for Fig. 2(a) can be obtained directly as follows, 



V 2 (a)(«) 



768ira q N r S 



5(iV c 2 -l)( Pl + P2 )4 « 



(24) 



where 



V. 



Ann 



99 



Aiia q 



(pi + P2) 2 
1 



(Ai - AD 2 



1 



6m 2 



1 



8iV c 

(9 " + 9 
^2 , 



(i - n • f 2 ) 



(^1 + <?2 



^2 



<?1 • ^2(9 + 



(25) 



where, q and q ' are relative momenta between quarks and antiquarks in the initial and final 



states, respectively. The isospin factor (1 — t\ ■ t^)/2 in (25) is introduced by considering 
the fact that the potential has a nonvanishing value only for isospin T = state of a quark 
and antiquark pair, which corresponds to the gluon quantum number. Performing Fourier 
transformation to V2( a )(9) yields 



48ira s N c S Ann 



5(iV 2 - l)m 



(26) 



where C/^ nn (x), the perturbative qq pair-annihilation potential in the coordinate representa- 



tion, is, 



5(t) 



7T 



(A 1 -AT) 2 (l-f 1 -f 2 ) 



4 16iV c m 2 
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_i_V')^)-A[ (?1 .v)(a 2 .V) 



-at ■ <t 2 V 2 



6(x] 



(27) 



The nonperturbative contributions of Fig. 2(b)-(d) may also be obtained by means of the 
same procedure as above. The total (/(/-pair annihilation potential can be obtained by sum- 
ming up the contributions of all diagrams including nonperturbative ones in Fig. 2 and the 
corresponding perturbative one, 



[4 nn(Total) (x) 



< nn (*) + 4 nn(NP) ( 



(28) 



where 



U. 



Ann(NP) 



qq 



(x) 



ira 



A8N C S | 1 m/(8m 2 - m? f )(0\q f q f \0) 
m 2 } 5(JV| - l)m 2 + W c y (4m 2 - m 2 ) 2 



Therefore, the effective potential between a quark and an antiquark of the same flavor may 
be expressed as the summation of the direct and annihilation potentials: 



U, 



qq of the same 



fl avo r (x)=^ CCt (^)+C n(TOtal) (^) 



(30) 



4. Nonperturbative QCD Corrections to the qq Pair-Excitation Potential 

Now, we turn to the nonperturbative QCD correction to a qq pair- excitation potential. 
By means of the same procedure as above, we can calculate the contribution of Fig. 3 to the 
effective potential. The result for Fig. 3(a) reads 



^><«> = w^w^M- (31) 

where V q ^ qqq ~(q), the usual perturbative quark excitation potential in the OGE approxima- 



tion, is [14 1 



V q ^ qqq {q) 



47ra s -^ 



l/i l 

2 \mi ni2, 



-q ■ (<7i x a 2 ) + 



, , , (32) 
with q = pi —pi. As suggested in Ref. [14|, we adopt two different approximations, i.e., for 



, Uq 



(case A), and uj q = 2rri2 with q ~ (case B). From the Fourier trans- 



formation, the expression for the transition potential (32) in the coordinate representation 
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can be written as 



and 



£j(B)q-»qqq^ 



i((Ti X <7 2 ) 



mi 



1 1 
+ | cr 2 

mi m 2/ 

x 2i&2 ■ V\ 



\x\ 2 m\ 



i6(t) AJAg 
2mJ 4 



Q~2 _^2_ 

mi m 2 



i((?l x <7 2 ) 



mi 



mi 



(33) 



(34) 



in case A and case B, respectively. In case A, the expression for ( |3l| ) in the coordinate 
representation is 



A? AS 



UrUx) = 6(t)-^4nai[D 3 \x\ d + D 2 \x\ 



(35) 



with 



Da 



8N r S 



5(N% - l)mi 



0"2 - Pi, 



(36) 



r> 2 



12iV,.5 



5(N* - 1 
In case B, one can easily obtain 



n • (ai x ga) 1 1 

h i{ 1 )(cr 2 ■ n) 

mi mi m 2 



48vra s AT r 5 



jy(B)q-+qqq ( 



X . 



(37) 



(38) 



Similarly, the potential for Fig. 3(b) does not give any contribution, i.e., U^t^ \x) = 0. 
The potential for Fig. 3(c) and that for Fig. 3(d) are the same, and turn out to be 



in case A, with 



U$(x) = U^ ) (x) 



(A) 



\a \a 



Ft\x\ +F + Y,Fo f)e ~ mfl " 
f 



(39) 



Pi '#2 (0|g/g/l0) 
AN c mi y rrif 



(40) 
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F n = 



8N r 



i ^<0|?/g/|0) 



m f 



mi 



-n ■ ((7l x «r 2 ) 





n ■ a 2 ) 


(-+-)< 


\mi 7712/ 





(41) 



and 



In case B, 



:-,(/) _ Pi -Q2 (0|g/g/|0) 
8iV c mi m) 

i (0|g/g/|0) 

16iV r 



1 1 \ / -* 

1 )(n ■ CT 2 

mi fni 



— n ■ (a i x cj 2 ) 
mi 



(42) 



E 



m/(0|g/g/|0) 



(4m2 - mj) 



1 + 



2(4m| 



m. 



(43) 



Therefore, the total transition potential in case A (or B) is 



Ur 



(A/B)q^qqq 



Total 



( X ) = ^(A/B) q ^qqq (:r) + ^(A/BJq-qqqCNP) ^ ^ 



(44) 



f/(A/B)q-»qqq(NP) non p er t U rbative correction to the perturbative qq-pair excitation 

potential from all diagrams shown in Fig. 3 is 



[/ (A)q^qq q (NP) ( 



in case A, or 

Lr (B)q^qqq(NP)/ 



X) 



+ U 3(c) ^ U 3(d) 



= 5(t)^^ s [D 3 \ 



xl 3 



+ D 2 \x\ 2 + 2Fi\x\ 



+ 2F + 2^F (/) e- m ^ |£| 
/ 



MB) + ^(B) + rfB) 



'3(a) 



y 3(b) "•" u 3(c) 

48N C S 



+ U. 



(B) 



3(d) 



5(iV c 2 - l)m 2 + V, 



1 m/(0jg/g/j0) 



(4m 2 . - 



ml 



1 + 



m. 



2(4m| 



(45) 



x[ /(B)q-qqq( a ,) ) 



(46) 
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in case B. 



5. Some Preliminary Numerical Results 

In order to provide new insight into the effect of the condensate correction on the quark- 
antiquark potential, we expect the obtained potential including the nonperturbative effect 
can be employed to improve the hadronic spectra and hadronic properties of J/^> and T 
families, for example the spin splitting between l 3 £i and l 1 ^ could be one of the sensitive 
quantities for the correction. Here we give the result for J/ty family, the nonperturbative 
effect in T family is not presented since it is qualitatively similar to that of J/^> family. 

Similar to Ref. ||, the Cornell potential U = — 4a s /3r + nr with corresponding 
parameters a s and k |l5|], which gave the best fit to the J/^ and T family data, is adopted 
as a basic condition, and the values of vacuum condensates are taken from Ref. ||. The 
total potential together with the phenomenological linear confinement can be written in the 
following form: 

U{r) = U C + Uf orT {r) + U? orr (r) + C/ 3 Corr (r), 
where JJi° rr is the correction from the non-trivial physical vacuum condensates, while 17% 1 



'2 

and U^ orr (sec Ref. ||) are the Breit-Fermi corrections to the Coulomb and linear confine- 
ment terms, respectively. 

The numerical calculation is performed in the following way. The Cornell potential is con- 
sidered as the dominant part of the potential, and the corresponding parameters a s and k, are 
determined before adding in the corrections. The values of a s and k are 0.381 and 0.182GeV 2 , 



respectively [15|. And then, the newly derived corrections due to the non- vanishing vacuum 
condensates as well as the Breit-Fermi corrections are treated as a perturbation adding onto 
the dominant part. The resultant values for the cc system are tabulated in Table 1. 
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Table 1 cc system 





exp't. 


U c = 

+ KT 

3r 


U c + U§ orr 


U c + Ui? orr 

+ U Corr 

P = 0.6 


U v + U^ orr 

JjCorr + jjCorr 
fjjCorr in Re f.g) 

P = 0.6 


U u + U^ orr 

JjCorr _|_ jjCorr 

P = 0.6 


I 1 s 


2978.8±1.9 


3074.0 


3010.8 


3046.9 


2979.3 


2981.3 


2 x a 


3594.0±5.0 


3662.1 


3625.0 


3646.4 


3446.3 


3466.5 


lV 


3096.88±0.04 


3074.0 


3095.1 


3083.0 


3090.6 


3094.1 


2 a si 


3686.00±0.09 


3662.1 


3674.5 


3667.3 


3493.1 


3516.1 


l a Po 


3415.1±1.0 




3440.8 


3410.5 


3321.2 


3330.7 


lV 


3510.53±0.12 


(1P C )3497.1 


3489.3 


3481.7 


3395.1 


3405.4 


1"P2 


3556.17±0.13 




3514.4 


3531.8 


3445.1 


3458.2 


E20 


141.07 





73.6 


121.34 


123.9 


127.5 


E21 


45.64 





25.1 


50.1 


50.0 


52.7 


aw 


118.08 





84.3 


36.1 


111.3 


112.8 




92.0 





49.5 


20.9 


46.8 


49.6 



In this table, £ 2 o = M l3p2 -M l3po , E 2 i = Mjs^-M^^, a£ = M^^-Mji^ and = M 2 3 Sl -M 2 i ao . 
(3, which characterizes the fraction of the confinement potential which comes from vector 
exchange, is taken as 0.6 ||. 

The experimental data are taken from "Partical Physics Booklet", July 1994, Partical Data Group. 



6. Discussion and Summary 

To have a deeper understanding of the hadronic structure, there have been various ways 
to modify the potential. As shown by Gupta et al. Q, Fulcher [16| and Pantaleone et al. 



|17], the nonperturbative condensate and the perturbative closed-loop corrections are two 
comparable effects, because these two kinds of modifications are in the same order of ct s . 
Hence, in a complete analysis, nonperturbative condensate corrections should be taken into 
account. 

The uncertainty due to the violation of the translational invariance in the fixed-point 
gauge is unavoidable. So as to include the nonperturbative corrections in a more self- 
consistent way in the nonrelativistic potential reduction, we used the covariant two-gluon 
VEV to evaluate the nonperturbative correction to the OGE potentials. 

In the resultant potentials, some new features appear: A linear term, which can play a 
role of the regular confinement potential, results from the condensate corrections of quarks 
and gluons to the gluon propagator; the modification of the gluon condensate to the gluon 
propagator can offer a |x| 3 term; The Yukawa-type term \x\~ 1 e~ mf in comes from the 
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nonzero quark condensate and may somewhat provide an interaction at longer range, such 
as part of the effects brought about by pseudoscalar meson exchange; the nonzero quark and 
gluon condensates can lead to the modification to the spin-orbital and tensor interactions, 
which may affect the baryon spectrum and scattering phase shift. It is noteworthy that the 
nonperturbative correction can not be expected to include as much as a purely phenomeno- 
logical ansatz, because this framework is an extrapolation from perturbative region. However, 
it does shed light on the physical picture and enrich our understanding of hadronic structure 
and the underlying mechanisms which determine how quarks are bound into hadrons. 

As a preliminary numerical result, the hadronic spectra of J/ty family and the spin 
splitting for this system were calculated by using the modified quark-antiquark potential, 
which indicates that the newly introduced corrections does lead to some improvements in 
the description of those hadronic observables to some degree. In the present work, only the 
lowest-dimension quark and gluon condensate corrections are taken into account, thus our 
calculation is reliable up to the intermediate range and it can only provide a description of 
the interaction between the quark and antiquark at the distance of r(r ~ 1 fm). For a longer 
distance scale, higher order condensates should be taken into account. 

In summary, it is the goal of this paper to extend the OGE potentials for the pure per- 
turbative region to the intermediate range in QCD and to investigate the gauge dependence 
of the QCD vacuum gluon fields. Here we only give an example to calculate them within the 
lowest condensate approximation and the result seems to be better than those obtained in 
the fixed-point gauge of the QCD vacuum gluon fields. In order to improve the numerical re- 
sult and give a reliable result for the higher excited states, the higher-dimension condensates 
should be included. Physically relevant results, such as the effective quark-quark interac- 
tion potential, should be gauge- independent. The differences between the present result of 
quark-quark interaction and that of Ref. H indicate that there are still gauge issues to be 
clarified. For this reason, we should fulfill the construction of description in Lorentz gauge 
by including the ghost propagators and ghost condensates as well as the (B 2 ) condensates 
according to the Slavnov- Taylor identities (STI) pjjfl . Further studies, such as the intended 
applications of the modified OGE potentials mentioned in the introduction and the improve- 
ment of description in Lorentz gauge, are in progress. Along this direction, we hope that 
the translational invariance and the gauge issues in the nonperturbative calculation can be 
clarified simultaneously. 
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Figure Captions 



Fig. 1. The Feynman diagrams for the contributions of the nonperturbative corrections to 
perturbative quark-quark potential in OGE approximation with the lowest dimensional quark 
and gluon condensates. 

Fig. 2. The Feynman diagrams for the contributions of the nonperturbative corrections to 
perturbative (/(/-pair annihilation potential in OGE approximation with the lowest dimen- 
sional quark and gluon condensates. 

Fig. 3. The Feynman diagrams for the contributions of the nonperturbative corrections to 
perturbative (/(/-pair excitation potential in OGE approximation with the lowest dimensional 
quark and gluon condensates. 
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